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Abstract 

Heavy right-handed neutrinos are of current interest. The interactions and decay of such neu- 
trinos determine their decoupling epoch during the evolution of the universe. This in turn affects 
various observable features like the energy density, nucleosynthesis, CMBR spectrum, galaxy for- 
mation, and baryogenesis. Here, we consider reduction of right-handed electron-type Majorana 
neutrinos, in the left-right symmetric model, by the WrWb channel and the channel originating 
from an anomaly, involving the SU(2)r gauge group, as well as decay of such neutrinos. We study 
the reduction of these neutrinos for different ranges of left-right model parameters, and find that, 
if the neutrino mass exceeds the right-handed gauge boson mass, then the neutrinos never decou- 
ple for realistic values of the parameters, but, rather, decay in equilibrium. Because there is no 
out-of-cquilibrium decay, no mass bounds can be set for the neutrinos. 



PACS numbers: 14.60.St,13.35.Hb,11.10.Wx,12.60.-i. 



1 E-mail address: paro27@hotmail.com 

2 E-mail address: allrc@usa.net 

3 E-mail address: amitava@cubmb.ernet.in 



1 



1 Introduction 



If a see-saw mechanism [pj is to account for left-handed electron neutrino v masses suffi- 
ciently small to be consistent with current ideas on neutrino oscillations, right-handed N 
neutrinos with mass, M, in the TeV scale come into the picture [|. High N masses in 
the range 1-20 TeV, and even higher masses, have been considered in studies of leptogenesis 
- baryogenesis [§, ||, || and e~e~ collisions Q. While some studies Q consider M smaller 
than M\y, the right-handed Wr boson mass, others specifically use N masses greater than 
the Wr mass p, p|. As cosmological and laboratory lower bounds for Z', Wr masses are 
of the order of 0.5 TeV, 1.6-3.2 TeV, respectively Q, there is no reason not to consider N 
masses greater than Z', Wr masses. Cosmological mass bounds for M, with M 3> Mw, 
have been considered in ||]-|l2|]. 

In a recent work |12| , (B + L)-violation from an anomaly, involving the SU(2)r gauge 
group |2|, |l3], p^l , was considered as a generation/reduction channel for N neutrinos satis- 
fying M > M\y. It was found that this anomalous channel played a role, in the decoupling 
of such neutrinos, at least as important as the NN — > channel (each of these chan- 

nels was found more important than the NN — > FF channel, F representing a relevant 
fermion). Matrix elements for NN — > FF and NN — > Wr were calculated in [Q from 
the Left-Right symmetric extension [15, 16| of the standard model, as an illustration. 

In the above work [12|, the N neutrinos were assumed to be stable, for simplicity. 
If, however, the Left-Right symmetric model (L-R model, hereafter) is to be taken as a 
serious working basis, N neutrinos cannot be considered to be stable. Decays involving 
v — N mixing, Wl — Wr mixing, and generation mixing, and CP-violating decays have 
been extensively studied. The last scenario has been widely used in generating lepton and, 
hence, baryon number from decay of Majorana-type neutrinos [Q, ||, 17, [l8|, [l9| . If the 
N mass is considered to be greater than the Wr mass, then these decay channels will be 
marginalised by the channel N — > Wr + e~ . Such a fast decay should have important 
effects on decoupling. It is this which is studied in detail in the present paper. The effect 
on leptogenesis has already received considerable attention p0|| . 

The effect of the decay of massive neutrinos on energy density, nucleosynthesis, the 
Cosmic Microwave Background Radiation (CMBR), galaxy formation, and stellar evolution 
has been well-studied |2l], In these studies, the decay time was taken to be large 

(typically, larger than 100-200 seconds), and the effect of decay at the crucial epochs followed 
up. The decay time was chosen greater than the decoupling time, i.e., it was assumed that 
the neutrinos first decoupled, and, then, their cosmological and astrophysical effects were 
felt, as a result of subsequent out-of-equilibrium decay [22[. If a fast decay like N — > W^ +e~ 
is considered, the relationship of decoupling and decay may not be like this, and the two 
have been considered, in the present paper, together, in one Boltzmann equation. 

The plan of the paper is as follows. Section 1 is the Introduction. In section 2, the 
thermally averaged NN — > W^W^ annihilation cross-section times relative velocity and 
the thermally averaged decay rate are calculated in the L-R model, assuming pure Majo- 
rana neutrinos. The anomalous rate per unit volume, used in ref.[12], is slightly modified 
to accommodate pure Majorana neutrinos. In section 3, the Boltzmann equation is writ- 
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ten down and treated approximately to obtain a decoupling criterion. The possibility of 
decoupling is investigated numerically. Section 4 discusses the conclusions. 



2 Annihilation, Decay, and Anomalous Reduction of Right- 
handed neutrinos: Thermal averages in L-R model 

2.1 Summary of relevant features of the L-R model 

In the L-R model, there are two doublets (v, ex) and (2V, exj) belonging to the representations 
(^,0, —1) and (0, |, —1), respectively, of SU(2)l x SU(2)r x U(1)b-l, where the quantum 
numbers refer to the values of Tl,Tr, and B — L respectively. To simplify the issues, only 
one generation is considered (the lightest), and N — u, Wl — Wr, and Z — Z' mixings are 
neglected. 

The symmetry-breaking from SU(2) L x SU(2) R x U(1) B -l -> SU(2) L x U(1) Y is 
achieved by means of a scalar triplet (T L = 0, T R = 1 and B - L = 2) A = (A++, A+, A ), 
putting < A >= vr/t/2 (chosen real). The right-handed gauge boson becomes massive 
due to the piece of the Lagrangian 

Law+w- = D^A^D^A (1) 

with = d^ + igT.Wfj, + i(g' /2)B^. From now on, the subscript R is dropped except when 
essential. The Tj form a 3 x 3 representation of the SU (2) generators in a spherical basis. 
Z' has the definition 

Z'^ = (Vcos 20/cos 6) W% - tan 8B^ (2) 
([]]) gives a W mass Myy = ^g 2 v R , and an interaction 

Lhw+w- = V2gM w W;W^H (3) 

where H(x)/V2 = Re(A°(x)) — vn/y/2, the displaced neutral Higgs field. 
Neglecting Z — Z' mixing, one gets the Z' neutral current piece 

L% f) = (g/2Vco S 28) cos ON^PrNZ'^. (4) 

6 is the weak mixing angle. The required charged current piece of the Lagrangian is 

L% e = (g/V2)(N^P R eW^ + h.c). (5) 

The Z' ,W + ,W~ interaction piece is 

L^ +w . = -ig(Vcos2e/ cos 0)[Z'HdpW+W v - - d u W+W u ~ - d^W~W v+ 

+d v w~w v+ ) + x -{d^z' v - d u zl)(w»+w- - w»~w v+ )] (6) 

The Majorana mass of N is thought to arise from the piece 

l rra = fR T Ce(T/V2).AR + h.c. (7) 
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where R T is the SU(2)r doublet (N, or), C is the charge conjugation matrix, and e = iT2- 
(r/\/2).A is the matrix 

( A+/^2 A++ \ 
I A -A+/^2 J ' 



and e(f/^/2).A is symmetric. This term gives a Majorana mass for the N neutrinos, M 
fvn/y/2, so that gj f = My/ /M = ryy, and an interaction 



L 



NNH 



2r w )N c HN 



(8) 



N c is the field conjugate to N. It has been pointed out that the effective one- 

loop mass matrix for multi-generation unstable Majorana neutrinos is not Hermitian and, 
strictly speaking, N 7^ N . But, here, at the level of tree-order calculations for a single 
generation, N = N c will be assumed. 



2.2 < a\v\ > for NN 



W R W R 



Three amplitudes %Mz> ,iM e ,iMn have been considered for NN — ► W + W , these arising, 
respectively, from NN — ► W + W~ in the s-channel, iViV — ► W + W~ in the t-channel, and 
NN —* W + W~ in the s-channel. Earlier calculations, with rjy <C 1 H, considered a heavy 
charged lepton exchange in the t-channel. In this paper, as required by the L-R model, an 
ordinary electron is considered to be exchanged. 

iMz> is calculated from (||) and iM e from (||), and iMn from (||) and (§). One finds 
that, in the limit M — > and s 3> M^-, iM e and iMz> cancel in tree-order, as expected. 
With a massive M 3> My/, an extra M-dependent term remains in iMz 1 + iM e . On 
calculating X) l-M| 2 > it is found that the interference terms between iMjj, and iMz' + iM e 
cancel, and for r\y <C 1, in the CM frame, 



w 



1 o n n 16 

¥ Eksm9cM+ M^ 



M 2 H ? 



0) 



(E, k) is the 4-momentum of the N neutrino in the CM frame, and 9qm, the angle of 
scattering, k has been written for |k|. 

The thermally averaged cross-section times relative velocity is 



< a\v\ > = (l/n eq ) J dir N dTTf l dTT W +dTT W -(2TT )5 (pn + Pn - Pw+ ~ Pw~) 
xJ2\M\ 2 e~ E »/ T e - E */ T . 
n eq is the equilibrium value of the number density n of the N neutrinos. 



n 



eq = 9N 



'MT 
W). 



_, 3 



-M/T 



(10) 



(11) 
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gN = 2 for Majorana neutrinos. The measure d-ni = gid 3 pi/{(2Ti) 3 2Ei]. |-A"f| 2 is the spin- 
averaged matrix element squared, with symmetry factor 1/2!, arising from the identification 
N = N c . 

The invariant integral / d7r w +dir w - (2tt) 4 5 4 (pn +Pn ~Pw+ ~ Pw~) 2~2 \-M\ 2 is calculated 
in the CM frame, then transformed to the comoving "lab" frame, in which N, N have 
energies En, E n , respectively, and thermally averaged according to (|l0|). In all calculations, 
N neutrinos are non-relativistic, viz. En = M + h N /(2M), k 2 ^ <C M 2 . This is because the 
interesting region for decoupling studies of a massive particle has T <C M. 

The result is 



< a\v\ >- 



1 ff 4 
2! 647rr^ 



r 16 

W + M 2 (4- 



,2 \2 



1 



3T(4 + r 2 H 
2M(4 - 




(12) 



r// = M H /M. 

The first term arises from Z'- and e-exchange, and apart from the 1/2! factor, agrees with 
the result of [1_1|, where the results of |J have been considered in the limit s — > AM 2 . 
The second term originates from //-exchange. In the calculation of this term, a further 



approximation has been made, 



viz. 



4k 2 < I4M" 



M||, 



z.e. 



this calculation is reliable 



provided Mjj is not very near in value to 2M. In Q, the if-exchange contribution was 
found to be negligible as s — » 4M 2 , because of a factor (s — AM 2 ) which arises for Dirac 
neutrinos. For Majorana neutrinos, this factor is absent, and this term cannot be neglected. 
As noted earlier, there is no interference between the //-exchange amplitude and those 



arising from Z'- and e-exchange. In [12], only the first term in ( |l2|) was considered, which 
is the approximation Mh 3> M. 

The processes NN — > FF, where F is a relevant fermion, have not been considered here. 
For M S> Mjy, the contribution of these processes is small compared to that of NN — > 
W + W~ |Q, and their effect on decoupling is overshadowed by the effects of iViV — ► W + W~ 
and the anomalous reduction of N neutrinos |Ol. 



2.3 Thermally averaged decay width of neutrinos 

The decay width for N — > W + + e~ can be calculated from (|5|). Calculation gives the 
spin-averaged matrix element squared in the neutrino rest frame 



with p = ±M(1 



W> 



1 ^P)MP 

2 P M 2 , 



being the momentum of the decay products, resulting in the width 



9 2 M 
32vr 



(1 



rlv? 1 + 



(13) 



In the frame in which the neutrino has energy En, the width becomes 



(M/E N )T e . 



(14) 
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The thermally averaged decay width [24, 25 1 is defined as 



r e = (l/n eq ) J dn N dn w+ dn w -(2n) 4 5\p N -p w+ - Pw -)J2\M\ 2 e~ EN/T 



= (l/n eq ) j (g N d 3 PN /(2irf)Tfe- E ^ T 
Using ( [Of ) and dlip, one gets 



(1 " f$r) 2 I 1 



32vr 



1 



2^. 



1 



3T 
2M 



(15) 



The calculation has been done in the approximation T <C M, <C M 2 . For M 3> M^, 



64vrr^ 



1 



3T\ 
2M J 



(16) 



If the mass of the N neutrino is greater than the A + mass, ([/]) allows the decay N — > 
A + + e~ . The corresponding 

J2\M\ 2 = (l/2)f(M 2 -Ml), 

where M + is the A + mass. In the rest frame of the neutrino, this decay width comes out 
as 

~2 icA-\r -\/r2 I%/t2 



T+ = (g z /647T)(M z /M^)M(l - r\ 
where r + = M + /M. There is another Yukawa piece of the Lagrangian 

L L Rt> = ^(hijLiQRj + h' ij L i T2^T 2 Rj) + h.c. 



(17) 



(18) 



The scalar field <3? transforms under the gauge group as (5,5,0) and is represented by the 
matrix 



■k'O 



Such pieces have been used in different models to break CP and induce baryogenesis through 
leptogenesis one generation is considered, the coupling constant h 

contributes to the electron mass and must be very small. To fit the observed baryon 
asymmetry, with N mass in the TeV range, \hij\ 2 values of the order of 10 -10 to 10~ 13 have 
been considered @, |5|. In this situation, as CP breaking has not been considered here, and 
only one generation taken into account, it has not been thought useful to consider iV-decays 
arising from (|l8|). In any case, it will be found that any enhancement of the decay width 
given in (|T5|) and (16) will only strengthen the main conclusion of the paper. 
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2.4 Anomalous generation of Majorana neutrinos 



The one generation, sphaleron-mediated, fermion number violating transition rate per unit 
volume, with |AL| = 1, \ AB\ = 1, for the quantum anomaly, involving the SU(2)r gauge 
group, was written, in [12], by extrapolation from the SU(2)l case [|(|, as 

21 7 /2 



Ar 



(1.4 x 10 6 ) 



x exp 




(19) 



M\y is the zero temperature Wr mass, z = Tr/Mw, where Tr is the critical temperature 
associated with the breaking of the SU (2)r gauge symmetry. So, z is essentially a quantity 
which reflects the uncertainty in the values of the L-R model parameters, while b captures 
the uncertainties involved in the extrapolation from the SU(2)l to the SU(2)r case in 



addition to those in the estimation of the prefactor of the anomaly driven transition [27|. 



In [12|, the anomalous rate of reduction of N neutrinos was considered, maintaining 
a distinction between N and N c . Here, N = N c neutrinos are considered, to maintain 
uniformity with the NN — ► W + W~ calculations. This entails an extra factor of 2, as seen 
below. 

For an anomalous process I, with Al = +1, such that 



l:i + j + 



N + a + b-\ 



one writes 12 



At 



dirNdTTad^b ■ ■ ■ dlTidlTj ■ ■ ■ \Mi\ 2 
X (2vr) 4 (5 4 (pAr +Pa+Pb~\ Pi - Pj )fNfafb 



(20) 



where I\ contains the result of the phase space integrations, apart from n eq [24, |25[, and 
i, j, ■ ■ ■ a,b, ■ ■ ■ are all supposed to be in equilibrium. 

Taking the view that leptogenesis and baryogenesis are effects of a smaller order, CP- 
symmetry is assumed. Then [O], for each process I, there is a AL = — 1 process 

V :N + a + b-\ ► i + j + ■ ■ ■ , 

with the same |.Mz| 2 . For this process, 



At, 



d-K^dlTadTib ■ ■ ■ dlTidTTj ■ ■ ■ 

x(27r) 4 5 4 {p N + Pa + p b + - 
I in 



Pi ~ Pj 



•)fNfafb ' 



(21) 
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The AL = +1 processes, generating N neutrinos, add up to 

Y, A i = n eq^k = -A R exp(-/3fj, L /2), 
i i 1 

from [26|. fj,i = [An is the chemical potential [p, e = 0, as electrons are in equilibrium). 
The 1/2 factor arises from the assumption [12], that, to a first approximation, the rate 



of generation of one member of a lepton doublet may be taken to be the same as that of 
the other, near equilibrium (not a bad condition at decoupling when the neutrinos are just 
falling out of equilibrium) . 
Similarly, 

= n Y. I i = \A R exp(+Pfi L /2). 
V i Z 



So, for small /i/y p6|. 



One gets 



n eq Ef k ~ ±A R (l-(3 m /2), 
nJ2iIi~ \A R {l+Pn N /2). 



^2h = A R /(n + n eg ), 
i 



and the anomalous rate of reduction of N neutrinos per unit volume 



An = J2 A P-E A i= ^-^Ar. (22) 



l n + n eq 



This has an extra factor of 2, compared to | |12[ |, because anti-particle processes, which had 
to be considered separately there, do not appear here, because of the assumption N = N c ' . 

3 Effect of decay on decoupling 

3.1 The Boltzmann equation for N neutrinos 

Using the results of the last section, one can write the Boltzmann equation 

(J i r j^ j 

— + 3Hn = -2 < a\v\ > (n 2 - n 2 ) - 2f e (n - n eq ) — (23) 

at H n + n eq 

where the second term on the left gives the effect of expansion, and the three terms on the 
right are to be taken from (0), (|l|), and (|22|), respectively. 



The expression for the Hubble parameter H is 

H = 1.66/ 5 T 2 /M Pl . (24) 

g* is taken ~ 100, and Mpi = 1.22 x 10 19 GeV. The 2 factor with < a\v\ > appears because 
two neutrinos are disappearing in NN — > W + W~ , considering N = N c ' . The 2 factor with 
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T e appears because of the two decay channels N — ► W + + e and N — > W + e .Writing 
x = M/T and y = n/s, where s = #* 5 (27T 2 /45)T 3 , with 100, (||) becomes 

y - y , 



ri)-d(i)(y-y e „) -<,(!) 



- eg 



y + y 



In (PD, 



/(*) 



1.41 x 10 ib GeV 



2.r 



(1 - rj./Af 



h(r H ) 



where a = l/rw = M/Mw and 

h(y) 



(i- 

d(x 
9(x) 



l + y 2 {(l-y 2 /8)(l-y 2 /4)+5/4} 



y 2 /4)3 

_ / 3.08 x 10 15 GeV N 
1 " ^ A% J 

( 3.06 x 10 23 GeV 



x exp 



M w 
-118.98 



ax 

)- 

/ a 



3 

1 

2x 



u- 



(25) 
(26) 

(27) 
(28) 



1 



ii- 



(29) 



with u = a/x = T/M< 



w ■ 



Writing Y = Y eq + A, it is noted that, before decoupling, Y is close to Y eq , and A' may 



be put equal to zero. Then, (25) may be put in the form 

A 



-Y' 

- 1 eq 



f( x )(2Y eq + A) + d(x) + 1 E 



From n eq in (|il|), one gets 



y, 



2.89 x lQ-^xae" 



• '</ • -• • and Y eq ~ 

at decoupling, when it is expected that x = 3> 1. The criterion for decoupling may be 
taken as A = c f Y eq , where d is of order 1. As in |p!2|| , c' is chosen to be 1. 
Then, the condition for decoupling is 

3f(x d )Y eq (x d ) + d{x d ) + = 1 (30) 

This is the key condition in our analysis of decoupling. Each factor on the lhs of (|3(i| ) is 
positive [28|. So, there is no scope for cancellation between different terms. We will show in 
the following subsections that, in fact, the condition can never be satisfied. The argument 
proceeds as follows. First, we consider (^), excluding the second (decay) term in the lhs. 
We show that there is a value x a , which we obtain numerically below, such that, for x < x a 
or x > x a , the lhs is, respectively, greater or less than unity in the absence of decay. We 
then check that for x > x a the decay term is much larger than unity so that there is no 
value of x (= M/T) for which eq. (30) is satisfied. 



-Y 

1 eqi 
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3.2 Decoupling in the absence of decay 



First, d(xd) is omitted, and the decoupling condition 

l(x) = 3f(x)Y eq (x) + = 1 (31) 

is solved to give x = x a . x a represents the value of M/T for which decoupling would occur 
in the absence of decay. 

Y eq gives a factor e~ x , and, so, the term with f(x) increases as x decreases. g(x) has 
an exponential factor of the form e~ Esp / T = e ~ Kx / a [ |l~2] , 
energy of the sphaleron mode which decays to cause anomalous L generation. The kinematic 
constraint on N production, E sp > M, gives K/a > 1 Jl^], and so g(x)/Y eq ~ e~(~~ 1 * )x , 
and, again, increases as x decreases. This means that l(x) > 1 if x < x a . So, x < x a will not 
satisfy the decoupling condition (^0|). Therefore, on the whole, it may be said approximately 
p4|| that, for x > x a , the annihilation rate plus the anomalous reduction rate is less than 
H. 

This programme is followed numerically. First, for definiteness, we fix b and z, the two 
parameters in the expression for g(x) - see eq. ( p9|) - at the values 1 and 4, respectively 
fllfl . a is varied from 2-100, for r#[] = 0, 1, 3, 10 and x a found by solving (|3l|) numerically. 



14, 26, 29, 30 



where E sp is the 



r H =l 


r H =3 


r H =W 


a 


X 


a 


X 


a 


X 


2 


25.79 


2 


25.06 


2 


22.30 


5 


28.49 


5 


27.74 


5 


24.89 


10 


30.54 


10 


29.77 


10 


26.86 


20 


32.60 


20 


31.81 


20 


28.83 


50 


35.32 


50 


34.52 


50 


32.18 


75 


73.77 


75 


73.77 


75 


73.77 


100 


238.83 


100 


238.83 


100 


238.83 



Table I: x = M/T at decoupling for different choices of a = M/Mw and rn — Mh/M without the 
inclusion of the decay contribution in the Boltzmann equation. Mw has been chosen to be 4000 GeV 



The results for Mw = 4000 GeV are shown in Table I for r#=l,3,10. There is a clear 
trend, showing an increase in x a for an increase in a (constant r#), and a decrease in x a 
for an increase in r# (constant a). The values for r#=100 and rjy=1000, we have checked, 
do not differ at all. It is safe to say that x a > 20 for the parameter ranges considered. 

We now address the uncertainty regarding the anomalous rate p7|]. The uncertainty 
is embodied in the parameters b and z appearing in the expression for this rate. The 
parameter z = Tr/Mw is of order unity (z = 3.8 in the SU{2)l case [^6|). To gauge the 

4 h(ru) assumes its asymptotic value of -2 for rn > 5. 
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a=2 


a=50 


6 


X 


b 


X 


10 y 


22.30 


10 3 


36.79 


10 2 


22.30 


10 2 


35.11 


10 


22.30 


10 


33.50 


1 


22.30 


1 


32.18 


lO" 1 


22.30 


io- 1 


31.60 


io~ 2 


22.30 


10~ 2 


31.47 


10~ 3 


22.30 


10~ 3 


31.46 



Table II: Effect of uncertainty in anomalous rate on the results of Table I 



sensitivity of the results on z, a calculation was made with a=50 (to give good weightage to 
the anomalous transition factor g(x)) and r#=l. x a , obtained by solving (pi), was 42 for 
z = 2, 35.6 for z = 3, and reached an asymptotic value of 35.3 for z > 4. For r#=10, x a was 
42 for z = 2, 35 for z = 3, and had an asymptotic value of 31.5. The conclusion x a > 20, 
of the last paragraph is not disturbed. In fll2]| , the uncertainty in b was taken into account 
by varying it through six orders of magnitude about 6 = 1. We do a similar calculation for 
z = 4, ru = 10, M\\r = 4000 GeV and present the results in Table II. It is clear that even 
this large variation of 6 through six orders of magnitude does not affect the conclusion of 
the previous paragraph. 

For comparison, the results for M\y = 2000 GeV, for the same values of the other 
parameters, are shown in Table III. 



r H =l 


r H =3 


r H =W 


a 


X 


a 


X 


a 


X 


2 


26.47 


2 


25.73 


2 


22.95 


5 


29.18 


5 


28.41 


5 


25.54 


10 


31.23 


10 


30.45 


10 


27.51 


20 


33.28 


20 


32.49 


20 


29.49 


50 


36.00 


50 


35.20 


50 


32.75 


75 


75.01 


75 


75.01 


75 


75.01 


100 


242.83 


100 


242.83 


100 


242.83 



Table III: Same as in Table I but for M w = 2000 GeV 



The expectation that l(x) - see eq. (|3l| ) - decreases below 1, as x increases through x a , 
was verified numerically for My/ = 4000 GeV, and a = 10,20,50,75, for each of r^=l,10. 
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Table IV shows the results for m = 1,10, and a = 10,50. 



r H =l 


r H =W 


a=10 


a=50 


o=10 


a=50 


X 


l(x) 


X 


l(x) 


a; 






l(x) 


20 


45.23 xlO 3 


30 


227.87 


20 


14.61 xlO 2 


25 


13.13 xlO 3 


35 


1.09 x 10~ 2 


40 


8.70 xlO~ 3 


30 


3.67 xlO" 2 


40 


1.51 xlO" 4 



Table IV: Behavior of l(x) ( |3l| ) around x a in the absence of decay. My/ = 4000 GeV for this Table. 



3.3 Effect of decay 

From the numerical results of the previous subsection it would be safe to say that for 
x < 20, l{x) > 1. Now, if one looks at d(x) (|28|), it is clear that, for x > 20, d(x) 3> 1, and, 
moreover, as x increases further, d(x) increases. So, there is no possibility of the decoupling 
condition (|30|) being satisfied. 



d(x), of course, has a simple physical meaning. In the lhs of (23), 



dn dY 
— — h 3Hn = Hsx—;—. 

dt dx 

Comparing (j23|) and (|25|), d(x) = (2F e /Hx). If d(x) > 1 for x = x a , with x a ~ 20, this 
means that F e 3> H at this point. Considering the physical meaning of x a , it may be 
concluded that, although, at temperatures lower than T a = M/x a , the annihilation rate 
plus anomalous reduction rate falls below the expansion rate, the decay rate remains much 
faster than the expansion, and this prevents decoupling. 

In the absence of decoupling, the fast decay constrains the N neutrino number density 
to follow the equilibrium density ~ e~ M l T . There is no out-of-equilibrium decay. 

One may check that this conclusion is not an artefact of the approximation M 2 ^> Myy 
in the calculations. Without this approximation, ( |i~5| ) gives the decay rate 



so that 



f = g 2 M(a 2 -l) 2 (a 2 + 2)(l-(3/2x)) 

64vra 4 1 ; 



„ . / 3.08 x 10 15 GeV \ (a 2 - l) 2 (a 2 + 2)x 
d(x) = — 



5 



2x 



(33) 



Simple calculation shows that if d(x) is to be < 1, for x ~ 20, then one must have a ~ 1, 
to 1 part in 10 7 (for M w = 4000 GeV). There is no reason for such fine tuning between N 
and W masses. 
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If the N decay width is augmented by the A + , e channel, i.e., if M > M + , 



d{x) 



'3.08 x 10 15 GeV 
M w 



(a 2 -l) 2 (a 2 + 2 



+ a(l-r 2 ) 2 







X 


V--} 




2s. 



(34) 



In this case, if d(x) is to be < 1 for x ~ 20, with My/ = 4000 GeV, not only must r\y ~ 1, 
but, also, r + ~ 1, each to 1 part in 10 7 , an unacceptable situation. 

Clearly, if the decay width is further augmented by introducing other channels, this will 
change T e , but not the x-dependence of f e or d(x), so that d(x) < 1 for x ~ 20 will be an 
even remoter possibility. 

d(x) can be decreased by increasing My/. Taking r# 
d(x) = 0.999 and l(x) = 0.001 for x = 4.25 and M w = 
x = 6.6 and M w = 1.6 x 10 17 GeV if a = 10. 



= 1, a little calculation shows that 
1.7 x 10 16 GeV if a = 2, and for 



This scale of My/ agrees with [2C] where the decoupling condition was, however, chosen 
simply as f e < H at T = M, f e having been assigned the value g 2 T /87T. In this paper, we 



have calculated T e 



327rM3 , and taken the thermal average of (M / E^)T e 

for T < M. We have used a decoupling criterion which includes the effect of annihilations 
and anomalous reduction, for a wide range of the parameters r# = Mh/M and r\y = 
My//M. Also, our method of approximate solution of the Boltzmann equation keeps the 
temperature of decoupling open and calculable, e.g., in the previous paragraph it was found 
that decoupling occurs for M > 10 16 , (10 17 ) GeV at T = M/4 (M/6), and not at T = M. 

So, the main point which emerges is that for there to be decoupling of massive neutrinos 
in the L-R model, it is necessary to consider right-handed gauge boson mass values far above 
the physically expected L-R mass scale. Further as M = aMyj , the see-saw mechanism 
will, then, give values of the v mass, which will be unacceptably small, when compared to 
neutrino oscillation values of Am 2 . 



4 Conclusions 

It is to be concluded that, in the L-R model, with right-handed neutrino mass greater 
than the Wr mass, even when the annihilation plus anomalous reduction rate for these 
neutrinos has fallen below the expansion rate, the decay remains faster than expansion, and 
becomes increasingly faster. So, massive right-handed neutrinos will decay while remaining 
in equilibrium. 

We find that this is true for a wide range of values of My/ /M and Mh/M, and also for 
a wide allowance of uncertainty in the anomalous rate. 

As the equilibrium number density varies as e~ M / T , right-handed neutrinos rapidly 
dwindle in number once the temperature falls below their mass. The decay products, with 
much lower masses, equilibriate immediately. Hence, there is no question of influencing the 
present density of the universe, CMBR, and nucleosynthesis or later events, and no mass 
bound can be set for right-handed neutrinos in the L-R model, if their mass is greater than 



the Wr boson mass (apart from an upper bound from unitarity jlC, |l 1| ) . 
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Because there is no out-of-equilibrium decay, the lepton (and baryon) number generation 
scenarios, utilising the decay of massive Majorana neutrinos will not work in the L-R model, 
if the neutrino mass is greater than the Wr boson mass, a result also previously noted in 
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